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The coupling with the lattice vibrations is shown to drastically modify the state-of-the-art picture 
of the excitonic states based on a frozen atom approximation. The zero-point vibrations renormalize 
the bare energies and optical strengths. Excitons acquire a non-radiative lifetime that decreases 
with increasing temperature. The optical brightness turns out to be strongly temperature dependent 
such as to induce bright to dark (and vice versa) transitions. The finite temperature experimental 
optical absorption spectra of bulk Si and hexagonal BN are successfully explained without using 
any external parameter. 
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The ab initio description of the excitonic states, ob- 
tained by solving the Bethe-Salpeter (BS) equation of 
Many-Body Perturbation Theory (MBPT), constitutes a 
well-established approach to interpret the photoexcited 
properties of bulk materials, surfaces, nanostructures and 
organic/bio-molecules Jj. Although absorption and pho- 
toluminescence (PL) experiments are usually performed 
at room temperature, in the standard approach the BS 
equation is solved assuming the atoms frozen in their 
crystallographic positions, thus neglecting the effect of 
lattice vibrations. As a consequence excitons turn out 
to be insensitive to the temperature T and to have an 
infinite lifetime. This is in stark contrast with the ex- 
perimental results, where the absorption and emission 
lines at any temperature show an intrinsic width that re- 
flects the finite lifetime of the underlying excitonic states. 
Moreover, in bulk semiconductors, it is a well known fact 
that the absorption line position, width, and intensity 
show a clear T dependence y. In the frozen-atom BS 
equation this dependence is not described at all. Even in 
the T ^ limit, where atoms vibrate to fulfill the un- 
certainty principle (zero-point vibrations) , the calculated 
absorption spectra is commonly convoluted with some ar- 
tificial, ad-hoc numerical broadening function chosen to 
yield the best agreement with the experiment. More gen- 
erally the finite temperature non-radiative damping, and 
the energy and optical strength renormalization, define 
the quantum efficiency of the excitons as photo-emitters, 
a key parameter in devising materials for optoelectronic 
applications. 

Bulk silicon (Si) and hexagonal boron nitride (/i-BN) 
are two paradigmatic semiconductors whose optical 
properties show remarkable differences. Si is one of 
the most deeply investigated material in the ab ini- 
tio community [l| . Although its optical properties have 
been studied in different theoretical frameworks [ij the 
finite temperature dielectric function, measured by Laut- 
enschlager and Jellison [§] twenty years ago, remains still 
unexplained. While Si has a small indirect gap, /i-BN 



is a wide direct gap quasi-two-dimensional semiconduc- 
tor [4]. At difference with Si, /i-BN optical spectra are 
dominated by a bound exciton with large binding en- 
ergy ('^O.TeV). Very recent results by Watanabe et al. Q 
stimulated interest in this material for the possible ap- 
plications as ultraviolet laser device. In the same exper- 
iment Watanabe observes a rich series of PL peaks that 
appear together with the free exciton bands pointing to 
a non-negligible effect of the lattice vibrations. 

In this Letter I solve, in a fully ab-initio manner, the 
Bcthe-Salpeter equation including the coupling with the 
lattice vibrations. The picture of the excitons obtained 
within a frozen-atom approximation turns out to be 
deeply modified, both at zero and finite temperature. Ex- 
citons acquire a non-radiative lifetime, otherwise infinite 
in the frozen-atom approximation. The finite tempera- 
ture optical spectra of Si and /i-BN are reproduced in 
excellent agreement with the experimental results. The 
thermal properties of the excitonic states are explained 
in terms of a weak (Si) and a strong (/i-BN) exciton- 
phonon coupling. In Si the lattice vibrations affect only 
the electron-hole substrate of the excitonic states, while 
in /i-BN, they participate actively in the excitons build- 
up. In /i-BN, this strong coupling induces bright to dark 
(and vice versa) transitions and reduces, at zero temper- 
ature, the lowest exciton binding energy by ~30%. 

In the frozen-atom (FA) BS equation the excitonic 
states |Afa} and energies are eigenstates and eigen- 
values of the Hamiltonian H^^, written in the elec- 
tron (e) hole (h) basis [1] 

i?f/ = {E, - En) Seh.e'h' + ife - h) ^ee' , (1) 
hh' hh' 

with -Ee(;i) and fe{h) the quasi-electron (hole) ener- 
gies and occupations. S is the Bethe-Salpeter kernel, 
that is a sum of a direct and an exchange electron- 
hole (e-h) scattering: E^^, = {eh\W - 2V\e'h'). W 

hh' 

is the statically screened and V is the bare Coulomb 
interaction. The absorption spectrum is given by the 



2 



imaginary part of the dielectric function, £2 {lo) 



where S^^ 



{GS\i£, ■ r\XFA) are the excitonic optical strengths, 77 
is a broadening parameter, V is the crystal volume, 
and ^ the light polarization direction. Eq.((T]) is ab 
initio because the single-particle quasi-energies and the 
kernel S are calculated starting from Density-Functional 
Theory (DFT) wavefunctions and energies, with no ad- 
justable parameters [bI]. In the standard approach, the 
quasiparticle (QP) energies i?e,/i, obtained within the 
GW approximation for the electronic self-energy [H , are 
assumed to be real and independent on T. This approxi- 
mation is justified by the fact that the smallest excitation 
energy in a semiconductor, the gap energy Eg, is usually 
much larger than the thermal energy corresponding to 
typical experimental temperatures, i.e. T <^ Eg/kb. The 
Hamiltonian H^"^ is then hermitian [ij and T indepen- 
dent. As a consequence, the energies Ef^^ are real, and 
rj is used as an a posteriori parameter to mimic the ex- 
perimental broadening of the absorption peaks . 

In the finite temperature regime the levels Ei acquire 
an explicit dependence on the temperature, Ei (T) — 
E, + AE, (T), with AE, (T) = AE^" (T) + ASf ^ (T). 
AE'^^ is the thermal expansion (TE) contribution [i,^. 
AE'^~P represents the complex energy correction that 
arises from the electron-phonon (e-p) interaction. In this 
work the e-p interaction is treated in the Heine, Allen, 
and Cardona approach Q where AE'^~p can be rewritten 
in terms of an e-p coupling function g^Fi (w). 



e—p 



(T) ^ JdLu g^E, (lu) [N {uj,T) + 1/2] , 



(2) 



with N (lo^T) — (^e^'^ — l) being the Bose occupation 
function. The complex g^Fi function is given by 



g'F,{u;)=J2 



dE, 



dN {lu,,T) 



5{uj - uj^) 



(3) 



with the sum is extended to all phonon modes v 
Eqs.(12HSl) tells that, 3? [AE'r^''] arises from the quadratic 
contribution to the expansion of Ei {T) — Ei in the atomic 
displacements. As shown in Ref. [11], the QP states can 
also decay emitting phonons, thus acquiring a finite life- 
time (oc 1/S [AiJj*^"^]). Phonons have very small ener- 
gies (« 100 me V) and, consequently, can be populated at 
the typical experimental temperatures. 

The temperature dependence of the QP states, aris- 
ing from the e-p interaction, modifies Eq.ljT]) 12]. As 
the Ef,^h {T) functions are complex, the BS Hamiltonian 
turns in a non hermitian operator 

Hee' (T) = i^f/ + [AE, (T) - AE,, (T)] 6,^,,,^' , (4) 

hh' hh' 

and the excitonic states are solution of the eigenprob- 
lem H{T)\\{T)) = Ex{T)\\{T)). The eigenstates 



IA(T)) are linear combinations of e-h pairs: |A(T)} = 
i:ehA^eh{T)\eh), with = {eh\\). If we plug this 
expansion in the definition of the excitonic energies 
i;A(T) = (A(T)]H]A(r)} weget 

ii;;,(r) = (A(r)|H^^]A(r))+ 

Y,\A^,^{T)\^[AE,{T)-AE^{T)]. (5) 

eh 

Using Eq.((5]) and neglecting the TE term, Eq.((S]) yields 

K [AEx (T)] = [(A (T) \H.^^\\{T)) - {\fa[R^^\\fa)] 
du K [g^Fx (w, T)] [N (uj, T) + 1/2] , (6) 



[Ex (T)] = / dec- 3 [g^Fx (a;, T)] [N (a;, T) + 1/2] , 



(7) 



where AEx (T) = Ex (T) - E^^ and I have introduced 
the exciton-phonon coupling function g'^Fx{u!,T) — 
Eeh\A^ehiT)\'[g^F,iu)-g^Ff,{Lo)]. Eqs.dMZD con- 
stitute a key result of this work. Eq.© defines, 
in an ab-initio manner, the non-radiative excitonic 
^ - [2Q=(£'a (r))]"\ that is otherwise 
FA approximation. The dielectric 
function now depends explicitly on T, e2(w,r) ~ 

and no damping 



lifetime 
infinite in the 



87T/Vj:^\SxiT)\^^[iu;~ExiT)y 
parameter r/ is needed anymore. 

Eq. ^ defines the temperature dependence of the exci- 
tonic energies and is composed of two contributions: the 
integral of the g'^Fx function arises from the renormal- 
ization of the electron-hole pairs \eh) that constitute the 
excitonic packet (with amplitudes A^fJ. This term rep- 
resents an incoherent contribution, where the electrons 
and holes interact separately with the lattice vibrations. 
The first term, instead, describes an active participation 
of the phonon modes in the excitonic state build-up. It 
is, then, a coherent contribution that modifies the 
components, and vanishes when \X{T)) = \Xfa)- Thus 
Eqs.([BH7]) define two physical regimes of the exciton- 
phonon interaction: in the weak coupling case ]A(r)) « 
\Xfa), and the incoherent contribution is dominant (this 
is the case of Si). /i-BN, instead, belongs to the strong 
coupling case, where the coherent term in Eq.® cannot 
be neglected. 

A remarkable property of Eqs.(IBH7]) is that although 
N {uj,T ^ 0) = 0, the excitonic energies do not reduce 
to the FA values and the excitonic lifetimes remain finite 
when T — > 0, because of the 1/2 factor. This factor arises 
from the quantum-mechanical vibrations of the atoms 
when T = (the so-called zero-point vibrations 13|). 

The experimental finite temperature optical spectra of 
Si [3], shown in Fig.(IT]), are dominated by two excitonic 
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peaks (El and E2) resonant with the electron-hole con- 
tinuum and characterized by a moderate e-h attraction. 
As the temperature increases, the Ei 2 peaks move to- 
wards lower energies, with a width that increases with 
T. This gradual red-shift has been studied only in an 
independent-particle approximation (IPA) [llj, thus ne- 
glecting excitonic effects. While the IPA shows only a 
weak dependence on T both peak position and widths of 
the Ei_2 peaks are well reproduced by the results of the 
finitc-T BS equation, shown in Fig.([T]). Excitons acquire 
a finite damping, that starting from ~ 30 meV at T = K 
and increasing to 60meV at room-temperature and 
~ 150 meV at T = 676 K, is in excellent agreement with 



the experimental estimations 11|. Compared with the 
frozen-atom BS equation the position of the Ei and E2 
peaks at T = is red-shifted by 80meV, so to correct 
the deviation of previous calculations Q from the exper- 
imental spectrum. 
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FIG. 1: (color on— line) Optical absorption of bulk Si (upper 
frames) for several temperatures and of /i—BN (lower frames) at 
room temperature. The experimental spectra [J Il5| (circles) are 
compared with the BS equation (solid line), and with the indepen- 
dent particle approximation (dot-dashed line). In the insets the 
exciton-phonon spectral functions 3? (^g^Fx {i^,T = 0)) are shown 
for the El (Si) and Bi (h-BN) peaks (see text). The width of the 
absorption peaks reflect the damping of the excitons due to the 
scattering with phonons. No additional numerical damping is in- 
cluded. The excitonic energies obtained within the frozen-atom BS 
equation (represented by the vertical dashed lines) are red-shifted 
in Si and blue— shifted in /i— BN. 



The g^F\ function can be now used to pin down the 



phonon modes that contribute to the red-shift of the Ei^2 
peaks. In the inset of Fig.((T]) the ^[g'^Fx (uj)) for the 
El state shows that the exciton is mainly coupled with 
the optical phonons (60meV peak), with the acoustic 
branches giving only a small correction. As the tempera- 
ture increases, the phonon population in Eqs.([5H7I) also 
increases, thus enhancing the red-shift and the width of 
the optical peaks, and leading to a linear scaling with 
the temperature when T > 200K and N {uj,T) - 1//3T. 
A more careful analysis of the different contributions to 
AEx (T) given by Eq.(I5]) shows that the incoherent con- 
tribution (second term in r.h.s. of Eq.(l6])) is dominant. 
This is due to the fact that the moderate e-h attrac- 
tion prevents the Ei_2 excitons to behave as a unique, 
bosonic-like, particles. Consequently the lattice vibra- 
tions mainly couple with the e-h substrate of the exci- 
tons. It is important to note that in this case, Eqs.(l6HZl) 
can be simplified using the result of the FA BS equa- 
tion, as AEx{T) w J dto [g^ Fx (uj)] [A (lj, T) + 1/2], 
with |A(r)) « \Xfa)- 
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FIG. 2: Temperature dependence of the energies and oscillator 
strengths of the near-gap excitons (bound and resonant) in /i-BN. 
The dashed line indicates the energy position of the optical gap in 
the independent— particle approximation. The size of the circles are 
proportional to the excitonic optical strength. The darker circles at 
T=0 is the result obtained neglecting exciton-phonon coupling. At 
T=0 we observe a 30% reduction of the Bi exciton binding energy 
due to the zero-point lattice vibrations. The B2 and _Ri excitons, 
instead, undergo a bright to dark (and vice versa) transition at 
room-temperature (see text). 

/i-BN is a anisotropic, insulating compound, consist- 
ing of graphite-like sheets with an hexagonal structure 
arranged in an ABAB. . . stacking The optical and 
electronic properties, as well as the lattice dynamics 14 1 
are strongly influenced by the layered structure. The 
in-plane experimental optical absorption spectrum mea- 
sured at room-temperature [l5| is shown in Fig. ([T]) , lower 
frames. Three prominent peaks are clearly distinguish- 
able: a bound state Bi at 5.98 eV, and two resonant 
states, Ri at 6.87 eV and R2 at 14.7 eV. The frozen- 
atom BS equation predicts the three peak energies to 
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be 5.75eV, 6.6eV, and 14.2eV [3, and 0.1 -0.5 eV red 



shifted if compared to the experiment. 

The room-temperature solution of the BS equation is 
compared with the experiment in Fig.([T]). Both experi- 
mental peak positions and widths are well described, and 
the Bi, Ri, and R2 states are blue-shifted of 0.07 eV, 
0.17eV, and 0.3eV compared to the frozen-atom BS 
equation results. The different sign of the phonon in- 
duced corrections of the excitonic peak positions is the 
first striking difference with the case of Si and can be 
understood looking at the function K {g^F\ (w)) for the 
Bi state, showed in the inset of Fig.((T]). The anisotropic 
structure of /i-BN is reflected in the rich series of phonon 
peaks in the g^F\ function. As showed by Serrano et 
al. [3], the phonon modes corresponding to the peaks 
at ~30meV and ~75meV are polarized perpendicularly 
to the hexagonal layers. As the bound excitons are spa- 
tially confined within the layer [l^, these modes tend to 
stretch the layers thus increasing the exciton localization, 
and consequently, its binding energy. The high-energy 
modes {uj > 100 meV), instead, are polarized parallel to 
the layer. These modes correspond to in-plane vibrations 
that interfere with the binding of the e-h pairs embodied 
in the excitonic state, counteracting the excitonic local- 
ization. Their stronger positive contribution to the g^F\ 
function causes an overall blue-shift of the absorption 
peaks, and a reduction of the exciton binding energy. 
Similarly to the case of Si, the /i-BN QP optical gap is 
shrank by the elect ron-phonon coupling by 0.12 eV. Thus 
we get an overall reduction of the lowest exciton binding 
energy of 0.2 eV, that is 30% the value obtained neglect- 
ing the exciton-phonon coupling (0.72 eV). 

The thermal evolution of the excitonic energies and 
optical strengths \Sx (T) ^ for the near-gap excitons is 
showed in Fig. . The size of the circles is proportional 
to \S\\'^. The opposite contribution to the g^F\ function 
of the low and high-energy phonons makes the excitonic 
energies to be almost constant for T < 500 K, in agree- 
ment with the experimental observation [l^. In contrast, 
the excitonic optical strength drastically depends on the 
temperature. We see that the Ri (resonant) and the 
B2 (bound) excitons undergo a bright to dark (and vice 
versa) transition at room-temperature. Indeed, we have 
that Sx (T) - = (G5'|e^«-^^[lA(r)) - \Xfa)], so that 
this astonishing effect is entirely due to the coherent con- 
tribution A [E\ (T)], given by the first term in Eq.®. A 
similar effects, due to the exciton spin decay, has been 
observed in quantum dots The present mechanism, 
instead, is independent on the dimensionality and driven 
by the temperature. 

From Fig. ([2]) we notice that the B2 and Ri transi- 
tions occur only when a bundle of states get close in en- 
ergy. The B2 state, for example, acquires optical strength 
only when it approaches the B3 state. The microscopical 
mechanism of the bright to dark (and vice versa) transi- 
tions is, then, a transfer of optical strength between en- 



ergetically close excitonic states. In the case of the B2 
exciton, for example, this process occurs by means of a 
mixing with the B3 state, that induces an increase of 
the contribution from bands with different parity in the 
e-h pairs embodied in the B2 state. This induces a fi- 
nite dipole, and a finite absorption cross-section. In the 
case of resonant excitons (with energy larger than the 
optical gap Eg) this hybridization is possible because of 
the continuum e-h substrate which connects the states. 
However, the bound excitons are discrete states, and the 
e-h substrate is replaced by the energy indetermination 
due to the finite damping. This is confirmed by a calcu- 
lation of the /i-BN excitons done neglecting the exciton 
damping, imposing the Ex (T) to be real. In this case the 
three bound states Bi^2,3 energies E\ (T) never cross, and 
the B2 state remains dark at all temperatures. 

In conclusion, the electron-phonon coupling induces a 
severe modification of the frozen-atom picture of the ex- 
citonic states both at zero and at finite temperatures. 
The proposed finite temperature Bethe-Salpeter equa- 
tion describes, in a fully ab initio manner, a wealth of 
new physical features, and makes clear that a proper and 
accurate description of the excitonic states in semicon- 
ductors and insulators cannot disregard the coupling with 
the lattice vibrations. 
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